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1. 


INTRODUCTION 


Similarity  of  the  velocity  profile  is  one  of  the  basic  concepts  of  fluid  flow  theory  going  back  to 
the  pioneering  work  of  Reynolds  [1]  in  the  1800’s.  Recently,  Weyburne  [2]  introduced  a  new 
approach  for  studying  similarity  of  the  velocity  profile  for  2-D  flow  along  a  wall.  The  following 
provides  a  more  concise  derivation  and  fixes  some  minor  omissions  from  that  report. 


Any  discussion  of  similarity  must  start  with  the  traditional  definition  of  similarity  of  the  velocity 
profile  given  formally  by  Schlighting  [3].  For  2-D  wall-bounded  flows,  the  definition  of 
similarity  is  that  two  velocity  profile  curves  from  different  stations  along  the  flow  are  similar  if 
they  differ  only  by  scaling  constants  in  y  and  u(x,y)  ,  where  y  is  the  normal  direction  to  the 

wall,  X  is  the  flow  direction,  and  L/(x,y)  is  the  velocity  parallel  to  the  wall  in  the  flow  direction. 


To  put  this  in  formal  mathematical  terms,  we  first  take  the  length  scaling  variable  as  and 
velocity  scaling  variable  as  .  These  scaling  variables  can  vary  with  the  flow  direction  (x- 
direction)  but  not  perpendicular  to  the  wall  (y-direction).  Hence,  according  to  Schlighting  [3], 
the  scaled  velocity  profile  at  a  station  along  the  wall  will  be  similar  to  the  scaled  profile  at  Xj 
if 


u(xi>yl^s{xi)) 


for  all  y. 


(1) 


From  a  mathematical  perspective,  it  is  found  that  this  definition  of  similarity  leads  to  a  number 
of  properties  of  the  scaled  profiles  that  can  be  used  to  advantage.  Consider  the  area  under  the 
scaled  velocity  profiles.  Similarity  necessarily  requires  that  if  the  scaled  profiles  are  plotted 
versus  the  scaled  y-coordinate,  then  the  area  under  the  curves  must  be  equal.  In  the  following, 
we  use  this  property  to  find  a  new  way  to  define  similarity  of  the  velocity  profile  as  well  as 
discover  some  properties  of  similar  velocity  profiles  of  2-D  wall-bounded  boundary  layer  flows. 

The  equal  area  concept  will  be  defined  in  terms  of  some  integral-based  requirements.  These 
results  are  derived  directly  from  the  definition  of  similar  curves  given  by  Eq.  1.  We  therefore  do 
not  need  to  know  how  the  curves  were  generated.  Hence,  the  flow  governing  equations  do  not 
have  to  be  invoked.  Furthermore,  the  new  results  apply  to  similar  profile  curves  whether  they 
are  laminar,  transitional,  or  turbulent  velocity  profiles  as  long  as  the  velocity  is  taken  as  the 
Reynolds  averaged  velocity. 
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2. 


SIMILARITY  OF  THE  VELOCITY  PROFILE 


2,1  The  First  Derivative  Profile 

We  start  by  considering  what  can  be  learned  from  this  equal  area  approach  to  similarity  for  2-D 
flow  along  a  wall.  In  the  analysis  below,  the  only  assumptions  that  are  necessary  as  to  the 
functional  form  of  the  velocity  profile  is  that  the  form  is  consistent  with  a  physically  realizable 
form,  i.e.  no  discontinuities  or  singularities.  The  development  presented  herein  is  based  on  a 
simple  concept;  for  similarity,  the  area  under  the  scaled  velocity  profile  curves  plotted  versus  the 
scaled  y-coordinate  must  be  equal  at  each  station  along  the  wall. 

Rather  than  starting  directly  with  the  consideration  of  the  velocity  profile,  we  will  first  consider 
the  implications  of  similarity  on  the  first  derivative  profile  curves  since  this  result  is  required 
later  on.  If  similarity  is  present  in  a  set  of  velocity  profiles  then  it  is  self-evident  that  the 
properly  scaled  first  derivative  profile  curves  (derivative  with  respect  to  the  scaled  y-coordinate) 
must  also  be  similar.  It  is  also  self-evident  that  the  area  under  the  scaled  first  derivative  profiles 
plotted  against  the  scaled  y-coordinate  must  be  equal  for  similarity. 


In  mathematical  terms,  the  area  under  the  scaled  first  derivative  profile  curve  is  expressed  by 


d[u[x,y/5^)lu^] 


where  o(x)  is  in  general  a  non-zero  numerical  constant.  For  clarity,  h(x),  nj(x),and 
have  been  shortened  to  h,  ,  and  .  Using  the  2-D  boundary  conditions  u(x,0)  =  0  and 
u[x,h)  =  L/g  (x)  where  (x)  is  the  free-stream  velocity  above  the  boundary  layer  edge,  and  a 
simple  variable  switch  (d {y  / =>  (l/ 5^  )dy) ,  Eq.  2  can  be  shown  to  reduce  to 


o(x) 

o(x) 

o(x) 

o(x) 

o(x) 


h/S, 


111 


d[u{x,y/5^)ju^] 


jdy 


d{u{x,y)/u. 


1  ^ 


dy 

du{x,y) 


S  0 


dy 


iF(x,y)] 


y=/7,0 


(3) 
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Similarity  at  two  stations  along  the  flow  located  at  xi  and  X2  requires  that 
o(x^)  =  o(x2)  =  constant.  Note  that  similarity  also  requires  that  h[x^)  j  ^  h[x2)  /  S^[x2) 

but  this  is  satisfied  as  long  as  h{x^)  or  h(x2)  are  chosen  to  both  satisfy  the  boundary  condition 
u{x,h)  =  Ug{x) .  Therefore,  for  similarity  of  the  velocity  profiles  for  a  2-D  boundary  layer  flow, 
the  velocity  scaling  factor  (x)  must  be  proportional  to  (x) ,  the  free-stream  velocity  above 

the  boundary  layer  edge.  Note  that  under  the  assumption  of  similarity,  Eq.  3  also  can  be  derived 
from  Eq.  1  directly.  The  above  derivation  is  offered  as  a  simple  introduction  to  the  equal  area 
approach. 


2,2  Defect  Profile 


Starting  with  the  formal  definition  of  similarity  given  by  Eq.  1  it  is  self-evident  that  for  the 
profiles  to  be  similar,  the  area  under  these  scaled  velocity  profiles  plotted  versus  the  scaled  y- 
coordinate  must  be  equal.  The  area  under  the  scaled  profile  curves,  in  integral  form,  is  given  by 


^o(^) 


1  UsJ 


(4) 


where  ho(x)  is  in  general  a  nonzero  numerical  constant.  Note  that  the  integral  is  written  using 
the  scaled  velocity  difference  rather  than  just  the  scaled  velocity.  It  is  simple  to  show  that 
similarity  of  the  defect  profile  u^-u  is  equivalent  to  similarity  of  the  velocity  profile  so  long  as 

Eq.  3  holds  true  (note  the  use  of  Eq.  3  was  not  explicit  but  subsumed  in  the  original  derivation 
given  in  Weybume  [2]).  The  use  of  the  defect  profile  has  two  advantages.  Eirst,  the  integral 
value  is  not  dependent  on  the  value  of  h  as  long  as  h  is  chosen  to  satisfy  the  boundary  condition 
u[x,h)  =  tVg  (x) .  This  simplifies  the  application  to  experimental  data.  Secondly,  using  the  same 
simple  variable  switch  from  above  and  simple  algebra,  Eq.  4  can  be  shown  to  reduce  to 


bo{x) 

bo{x) 

bo{x) 


h/Ss 

j.. 


y  \  u^-u[x,y/5^) 


1  ^ 


Ue-u{x,y) 


S  0 


bo{x)  = 


Us  i 

dlU, 


u^-u{x,y) 


S.  Uc 


where  the  is  the  displacement  thickness  given  by 


S,{x)  =  jdy  {l-u{x,y)/u^}  , 


(5) 


(6) 
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and  where  h  satisfies  the  boundary  eondition  u[x,h)  =  (x) .  Eq.  5  is  an  exaet  equation  that 

applies  whether  the  profiles  are  similar  or  not.  Similarity  at  two  stations  along  the  flow  loeated 
at  xl  and  x2  requires  that  bg  )  =  't’o  (•’^2 )  =  constant .  Note  that  it  is  also  neeessary  that 
h{x^)/S^{x^)  =  h{x2)/S^{x2)  and,  as  above,  either  h{x^)  or  h(x2)  can  be  freely  chosen  as 
long  as  they  both  satisfy  the  boundary  condition  u[x,h)  =  (x) .  The  importance  of  Eq.  5  in 

regards  to  similar  profiles  is  that  it  means  that  the  thickness  scaling  factor  and  the  velocity 
scaling  factor  cannot  be  independent  for  2-D  wall-bounded  similarity  flows. 

Combing  the  result  given  by  Eq.  5  with  the  result  given  by  Eq.  3,  then  it  is  evident  that 


c(x) 


S^{x) 


(7) 


where  c(x)  is  in  general  a  non-zero  numerical  constant.  Similarity  at  two  stations  along  the  flow 
located  at  xi  and  X2  requires  that  c(x;^ )  =  c(x2 ) .  Eq.  7  is  important  in  that  it  states  that  if 
similarity  exists,  then  the  displacement  thickness  must  be  a  length  scale  that  results  in  similarity. 


Having  equal  hg  (x)  values  at  different  stations  along  the  flow  is  a  necessary  but  not  a  sufficient 
condition  for  similarity  of  a  set  of  profile  curves.  If  the  scaled  velocity  profiles  are  similar,  then 
it  is  self-evident  that  the  scaled  velocity  profile  curves  multiplied  by  the  scaled  y-coordinate 
raised  to  the  nth  power  must  also  be  similar.  In  area  integral  terms,  the  area  under  the  scaled 
velocity  profiles  multiplied  by  the  scaled  y  -coordinate  raised  to  the  nth  power  is  equivalent  to 


bn{x) 


y 

1  UJ 

y 

u^-u{x,y/5,) 


(8) 


where  h„  are,  in  general,  non-zero  numerical  constants.  It  is  self-evident  that  a  sufficient 
condition  for  similarity  at  two  stations  along  the  flow  located  at  xi  and  X2  is  that 

h„(xi)  =  h„(x2)  forn=0,l,2,3,  ...,Go  ,  (9) 

so  long  as  h  is  chosen  appropriately  as  discussed  above. 
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3. 


DISCUSSION 


An  important  point  about  the  above  derivations  is  that  although  the  results  are  not  presented 
formally,  the  results  above  are  mathematieally  rigorous  and  ean  be  easily  substantiated  in  the 
form  of  mathematieal  proofs. 

The  approaeh  used  to  generate  these  new  results  is  different  than  the  traditional  method.  In  the 
traditional  method,  one  starts  with  the  flow  governing  equations,  then  one  makes  eertain 
assumptions  to  reduee  the  governing  equations  to  ordinary  differential  equations  and  tries  to  find 
appropriate  analytical  expressions  for  the  solutions.  Hence  the  traditional  method  is  intimately 
bound  to  the  physics  of  the  flow  through  the  flow  governing  equations.  There  may  be  a  tendency 
to  dismiss  the  new  approach  described  above  since  the  flow  governing  equations  are  not 
invoked.  However,  note  that  one  defines  similarity  of  a  set  of  velocity  profile  curves  using  Eq.  1 
[3].  This  equation  is  asking  a  simple  question;  are  the  left  and  right  sides  equivalent  or  are  they 
not.  There  are  no  flow  governing  equations  invoked  in  posing  this  question  or  answering  the 
question.  If  one  experimentally  measures  the  velocity  profile  at  two  stations  along  the  flow,  it  is 
not  necessary  to  invoke  the  flow  governing  equations  to  test  if  the  scaled  profile  curves  are 
similar.  That  is,  the  similarity  solutions  to  the  Navier-Stokes  equation  are  not  similar  because 
they  are  derived  from  the  Navier-Stokes  equation;  they  are  similar  if  one  substitutes  the  velocity 
solutions  from  the  Navier-Stokes  equation  into  Eq.  1  and  then  finds  that  the  scaled  velocity 
profile  curves  are  equivalent  at  different  stations  along  the  flow. 

The  similarity  scaling  result  given  by  Eq.  5  is  important  in  that  it  states  that  if  similarity  exists, 
then  the  similarity  length  scale  factor  and  the  similarity  velocity  scale  factor  cannot  be 
independent.  The  form  of  this  interdependence  is  new  and  different  from  than  that  found 
previously.  This  interdependence  is  not  like  pressure  gradient  parameter  A  found  in  Castillo 
and  George  [4]  for  example.  The  pressure  gradient  parameter  A  provides  a  required  relationship 
for  similarity  between  the  boundary  layer  thickness,  the  velocity  at  the  boundary  layer  edge,  and 
the  gradient  of  these  two  parameters  in  the  flow  direction.  The  new  interdependence  factor  given 
by  Eq.  5,  on  the  other  hand,  provides  a  necessary  similarity  relationship  on  the  proposed  scaling 
factors  in  terms  of  two  experimentally  accessible  variables.  However,  the  biggest  difference 
between  the  two  parameters  is  that  whereas  Eq.  5  is  invoked  directly  from  the  definition  of 
similarity  (Eq.  1),  the  Castillo  and  George  parameter  is  based  on  looking  at  the  x-momentum 
balance  equation  in  the  outer  region  of  the  turbulent  boundary  layer.  Note  that  if  one  looks  at  the 
inner  region  of  the  turbulent  boundary  layer,  a  different,  much  more  restrictive  parameter  is 
obtained.  Castillo  and  George,  nor  anyone  else,  has  explained  why  it  is  okay  to  ignore  the  inner 
region  restrictions  when  dealing  with  velocity  profile  similarity  using  their  approach.  In 
contrast,  there  are  no  restrictions  in  the  approach  developed  herein.  It  applies  to  wall-bounded 
turbulent  boundary  layers,  laminar  flows,  and  transition  flows.  If  similarity  exists  in  any  of  these 
flows  then  the  results  described  above  must  hold  true. 

One  remarkable  finding  is  that  if  we  take  in  Eq.  5  as  the  boundary  layer  thickness  S ,  then 
(x)  (with  hg  =  constant )  is  the  empirically  derived  velocity  scale  developed  by  Zagarola  and 
Smits  [5]  for  turbulent  boundary  layer  flow.  Zagarola  and  Smits  and  others  have  shown  that  the 
velocity  scaling  factor  given  by  Eq.  5  with  as  the  boundary  layer  thickness  can  collapse 
certain  experimental  turbulent  profiles  to  a  single  curve  (i.e.,  they  behave  similarly).  The  above 
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results  now  give  a  solid  theoretieal  foundation  to  the  empirical  results.  Note  that  although  this 
equation  explains  the  empirical  approach  of  Zagarola  and  Smits,  it  is  important  to  point  out  that 
the  fact  that  the  similarity  length  scale  factor  and  the  similarity  velocity  scale  factor  must  follow 
Eq.  5  (with  (x)  =  constant )  is  not  something  appreciated  in  the  fluid  flow  community.  Indeed, 
if  one  checks  the  recent  literature  dealing  with  similarity  of  wall-bounded  turbulent  flows  [4-6], 
it  is  apparent  that  no  one  has  made  the  connection  that  similarity  requires  that  the  length  and 
velocity  similarity  variables  must  be  coupled  through  Eq.  5. 

By  far  the  most  significant  new  finding  is  that  for  the  thickness  scaling  variable,  Eq.  7  indicates 
that  if  similarity  exists  in  a  set  of  velocity  profiles,  then  the  displacement  thickness  must  be  a 

similarity  length  scale  parameter.  This  new  result  applies  to  all  2-D  wall  bounded  boundary 
layer  flows.  This  is  the  first  time  that  a  particular  parameter  has  been  identified  as  a  similarity 
length  scale  and  most  significantly,  it  is  the  first  time  that  the  displacement  thickness  has  been 
identified  as  a  similarity  length  scale  parameter.  The  importance  of  this  result  cannot  be 
overstated.  Consider  that  the  2-D  wall-bounded  similar  flows  encompass  the  Blasius  [7]  and  the 
Ealkner-Skan  [8]  similarity  solutions.  In  the  more  than  100  years  since  Blasius  published  his 
work  no  one  has  ever  made  the  connection  that  the  displacement  thickness  is  a  similarity 

length  scale  for  the  Blassius  flow  solution.  The  same  is  true  for  the  Ealkner-Skan  similarity 
solution.  Up  until  now,  the  length  scale  for  the  Blasius  solution  has  been  taken  as  the  square  root 
of  the  kinematic  viscosity  times  the  distance  along  the  wall  divided  by  the  velocity  at  the 
boundary  edge  (infinity).  The  connection  to  the  displacement  thickness  has  never  been  made 
theoretically  but  was  posited  as  a  possibility  by  Weybume  [9]  based  on  examining  certain 
experimental  data  sets. 

Another  finding  from  above  is  that  Eq.  3  indicates  that  the  appropriate  velocity  scaling  constant 
for  2-D  wall  bounded  flow  similarity  is  .  We  note  that  Castillo  and  George  [4]  came  to  a 
similar  conclusion  using  a  momentum  balance  type  approach  to  similarity  for  flows  with  a 
pressure  gradient.  In  contrast,  the  results  above  apply  to  all  2-D  wall  bounded  flows  displaying 
similarity  regardless  of  whether  a  pressure  gradient  exists  or  not. 

It  must  be  emphasized  that  the  results  obtained  above  do  not  exclude  other  parameters  from 
being  similarity  scaling  variables.  Consider  for  example  the  Prandtl  Plus  scaling  variables. 

Jones,  et  al.,  [6]  used  a  flow  governing  scaling  approach  to  show  that  for  the  case  of  zero- 
pressure-gradient  turbulent  boundary  layer  in  the  limit  of  infinite  Reynolds  number,  the  flow 
friction  velocity  leads  to  a  valid  similarity  solution.  Eor  the  non-infinite  limit  case,  we  can 

postulate  that  is  a  possible  similarity  variable.  That  would  mean  that,  based  on  Eq.  3, 


^e(^) 

u,(x) 


constant 


(10) 
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This  is  of  course  the  well-known  Rotta  [10]  condition  stating  that  for  similarity,  the  inner  and 
outer  regions  veloeity  sealing  variables  must  evolve  at  the  same  rate  as  one  moves  along  the  wall 
in  the  flow  direetion.  What  is  interesting  here,  is  that  we  ean  now  also  establish  a  length  seale 
ratio  that  must  hold  if  the  Prandtl  sealing  variables  are  to  result  in  similarity,  that  is 


V 


constant 


(11) 


This  means  that  if  similarity  is  present  in  set  of  velocity  profiles  and  if  the  Prandtl  scaling 
variables  are  to  be  similarity  sealing  variables,  then  Eqs.  10  and  1 1  must  apply. 

Mathematically  it  is  evident  that  a  set  of  arbitrary  curves  are  similar  if  the  condition  given  by  Eq. 
9  exists  for  all  n.  However,  for  flow  similarity  of  veloeity  profiles  taken  at  various  stations  along 
the  wall  in  the  flow  direetion,  the  profile  curves  are  not  arbitrary  and  so  it  may  only  be  neeessary 
to  insure  that  ho  (•’^i )  =  (•’^2 )  •  case,  this  approach  to  similarity  has  an  advantage  from  an 

experimental  standpoint  since  the  equal  area  test  method  would  allow  for  statistieal  testing  for 
similarity  by  comparing  (x)  values  at  various  stations  along  the  wall.  Erom  a  praetieal 
standpoint,  this  method  is  superior  to  the  use  of  Eq.  1  sinee  to  use  Eq.  1  the  experimentalist 
needs  to  insure  the  veloeity  at  eaeh  y  / (x)  value  is  equal  at  eaeh  measurement  station.  This  is 

a  very  diffieult  task  sinee,  in  general,  the  aetual  values  of  (x)  are  not  known  a  priori.  As  a 
result  of  the  y  /  (x)  issue,  the  usual  imprecise  method  the  flow  community  presently  uses  to 

judge  whether  a  set  of  velocity  profiles  are  similar  is  to  plot  the  profiles  (with  lines  eonnecting 
the  data  points)  and  use  the  subjective  “chi-by-eye”  method  to  judge  the  suecess  or  failure.  If  the 
plotted  curves  fall  on  top  of  one  another,  they  are  considered  similar.  In  contrast,  the  equal  area 
test  method  would  allow  for  statistieal  testing  by  performing  simple  numerieal  integrals  (Eq.  8) 
of  the  profile  data  while  insuring  Eqs.  3  and  7  are  constant  at  the  stations  along  the  flow. 
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4. 


CONCLUSION 


A  new  approach  for  studying  similarity  of  the  velocity  profile  was  outlined.  It  starts  from  the 
equation  used  to  define  similarity  of  the  velocity  profde.  This  method  was  used  to  discover 
fundamentally  new  results  for  the  similarity  of  the  2-D  wall  bounded  boundary  layer  velocity 
profiles.  It  was  shown  that  if  similarity  exists,  then  the  similarity  velocity  and  length  scaling 
constants  cannot  be  independent.  Furthermore,  it  was  shown  that  if  similarity  exists,  the 
displacement  thickness  must  be  a  length  scaling  variable  and  the  velocity  at  the  boundary 
layer  edge  must  be  a  velocity  scaling  variable. 
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